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Abstract 

We study the phase space of spatially homogeneous and isotropic cosmology 
in general scalar-tensor theories. A reduction to a two-dimensional phase space is 
performed when possible — in these situations the phase space is usually a two- 
dimensional curved surface embedded in a three-dimensional space and composed 
of two sheets attached to each other, possibly with complicated topology. The 
results obtained are independent of the choice of the coupling function of the 
theory and, in certain situations, also of the potential. 
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1 Introduction 



Brans-Dicke and scalar-tensor gravity have been studied for over four decades. The 
original motivation for the introduction of Brans-Dicke theory [1] came from Mach's 
principle and cosmology but since Brans and Dicke's original paper, significant motiva- 
tion has been added for research in scalar-tensor gravity. First, Brans-Dicke theory has 
been generalized to scalar-tensor theories with coupling functions instead of coupling 
constants [2]. More important, it has been realized that scalar-tensor gravity contains 
features that are common to supergravity and string or M-theory. The low-energy limit 
of the bosonic string theory reduces to a Brans-Dicke theory with Brans-Dicke parameter 
u = -l [3]. 

Renewed interest in cosmology in scalar-tensor gravity (see Refs. [4, 5]) has been 
motivated by the extended [6] and hyperextended [7] scenarios of inflation in the early 
universe and by scalar-tensor models of dark energy trying to explain the present acceler- 
ation of the universe [8] . Even with the simplifying symmetry of spatially homogeneous 
and isotropic Friedmann-Lemaitre- Robertson- Walker (hereafter "FLRW") cosmology, 
the field equations of scalar-tensor gravity for the metric and the Brans-Dicke-like scalar 
field are complicated, only a few exact solutions are available, and a phase space pic- 
ture is very useful for the qualitative description of the dynamics. A few excellent and 
detailed works on the dynamics of Brans-Dicke homogeneous and isotropic cosmological 
models are available [9]— [11]. The possibilities studied in the literature include giving 
a mass or other potential V((f>) to the Brans-Dicke-like field (as is natural in particle 
physics), adding a perfect fluid while keeping V = 0, adding a cosmological constant A, 
and/or considering FLRW spaces with the three possible values 0, ±1 of the curvature 
index K. However, the geometric structure of the phase space, including its dimension- 
ality, is still poorly understood for scalar-tensor theories more general than Brans-Dicke 
theory. In addition, the studies of the dynamics available in the literature often use 
variables that mix the scale factor a(t) of the FLRW metric with the Brans-Dicke-like 
field 0. Although convenient from the formal point of view, the use of these variables 
tends to obscure the physical picture and it may be difficult to rephrase results in terms 
of the physical variables traditionally used in cosmology, such as the scale factor a, the 
Hubble parameter H, the comoving time t, and the scalar field 0. 

In this paper we focus on spatially homogeneous and isotropic FLRW spaces de- 
scribed by the metric 



ds 2 = -dt 2 + a 2 (t) 



<h ' 2 + r 2 (d9 2 + sin 2 9 d<p 2 ) 



(1.1) 



in comoving coordinates (t, r, 9, (p) and with curvature index K. We consider scalar- 
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tensor theories of gravity described by the action 




0i?-^^ b V a 0V fe 0-\/(0) 



+ S (m) , (1.2) 



where is the action describing ordinary matter, g is the determinant of the metric 
tensor g a b, V c is the covariant derivative operator associated with g a b, R is the Ricci 
curvature of spacetime, and (f> is the gravitational Brans-Dicke-like scalar self-interacting 
through the potential V ((/>). The metric signature is — , +, +,+,□ = g fec V&V c , and we 
follow the notations and conventions of Ref. [12]. 

The study of FLRW cosmology at very early times has limited value because isotropiza- 
tion is probably achieved by inflation, which washes away information about the pre- 
inflationary state of the universe. At very early times before inflation the universe may 
have been anisotropic (see Refs. [13] for studies of isotropization of Bianchi models in 
scalar-tensor gravity). However, it is meaningful to consider FLRW universes during and 
after inflation and in particular during the present and future evolution of the universe. 
In the quest for successful models of dark energy able to explain the present acceleration 
of the universe deduced from the study of type la supernovae [14] , scalar-tensor gravity 
has been used many times [8]. In this context, it is suggested [15] that the universe may 
evolve into a chaotic regime, thus making it impossible to predict its evolution in the 
distant future. The dimensionality of the phase space is important in this regard because 
a dimensional reduction of the phase space to two dimensions (which is achieved in this 
paper for many situations in the context of general scalar-tensor theories) is believed to 
inhibit the presence of chaos (see Ref. [16] for the special case of a nonminimally cou- 
pled scalar field theory). The impossibility of chaos in a two-dimensional phase space 
is guaranteed when the phase space is flat but it is not entirely trivial when the phase 
space is a curved two-dimensional surface, as is often the case in scalar-tensor cosmology 
[17]. 

In this paper we do not provide a complete phase space analysis as is done in previous 
works for Brans-Dicke theory and for special choices of the scalar field potential V(<t>) and 
of the coupling function w(<j>) [9]— [11], but we focus instead on the geometric structure 
and dimensionality of the phase space for any coupling function u(4>) and (except special 
cases) for any potential V((j>). A complete phase space analysis cannot obviously be 
performed without specifying the functions u and V, but we believe that a unified 
treatment of the phase space geometry in the general case is of interest, especially in 
view of the discussion of chaotic dynamics and of the uncertainty in the form of uo and 
V. Results similar to those obtained in this paper have been presented elsewhere for a 
special scalar-tensor theory [18]-[22], that of a scalar field coupled nonminimally with 
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the Ricci curvature and described by the action 



g(NMC) 



d 4 x 



-9 



+ s (m) , 



(1.3) 



where k = 8ttG and £ is a dimensionless coupling constant (in our notations £ = 1/6 
corresponds to conformal coupling). This theory is formally a scalar-tensor theory. In 
fact, by redefining the scalar field according to 



G 



the action (1.3) is written as 



g(NMC) 



I 



d x 



16tt 



ipR 



g ab V a ipV h ip-U{ci>) 



+ S (m) 



where 



and 



u(<p) 



4£ (1 - dp) 



U(<p) = 16ttV [<j>(ip)\ = 16nV d= 1 



G<p 



8ttG£ 



(1.4) 



;i.5) 



;i.6) 



(1.7) 



The theories in which the geometry and the dimensionality of the phase space of scalar- 
tensor FLRW cosmology are known are summarized in Table 1. This table summarizes 
also the main results of the present paper, which are derived in the following sections. 



2 The equations of scalar-tensor cosmology 



In this section we consider scalar-tensor gravity as described by the action (1.2) and we 
write the field equations in a form that is convenient for later use. The field equations 
obtained by varying the action (1.2) are 



G 



ab 



(f) 2 



V a 0V fe 0- -g ab V c <j)V c 



+ i (V a V b - g*U4>) - ^g ab + j Tff , 

(2.1) 



^=2^(^-^-^^ + 8 ^ M )' 



(2.2) 



3 



where is the stress-energy tensor of ordinary matter (usually treated as a fluid) 
which is covariantly conserved, V 6 T ( [™' ) = 0, and T = T a a . Note that the combination 
<f)dV/d<f) — 2V in eq. (2.2) disappears if the potential consists of a pure mass term V((f>) = 
m 2 <j) 2 /2 (which is the special case considered in Ref. [10] for Brans-Dicke theory). In the 
FLRW metric (1.1) the field equations assume the form 




H' = -H\r\+-^.IZ \ - + -£—, (2.3) 



H 




K 



8tt 



a 2 (2uj + 3) . 



[(u + 2) p^ + uoP^] 



^-2V + - 

d(j) d(j) 



(2.4) 



2u + 3 



2V-0^ + 8vr (p^ 



3P (m) ) 



(2.5) 



where an overdot denotes differentiation with respect to the comoving time t, H = a/a 
is the Hubble parameter, and p^ and p( m ) are the energy density and pressure of the 
material fluid, respectively. The latter satisfies the conservation equation 



p (m) + 3H (p (m) + P (m) ) = . 

We assume that the matter fluid has equation of state 

p{m) = ( 7 _ i) p H 

with 7 a constant. Then eq. (2.6) is immediately integrated to yield 

(m) _ _Po_ 



a 3 T 



(2.6) 



(2.7) 



(2.8) 



(with p a constant), a fact that is used in the following sections. In the general case, the 
dynamical variables are a and 0, the dynamical equations (2.3) and (2.5) are of second 
order in a and (there are only two independent equations in the set (2.3)-(2.5)), and 

the Hamiltonian constraint reduces the dimensionality of the phase space (a, a, 0, (j^j to 

three, as noted in Ref. [10] for the special case of Brans-Dicke theory. In many situations 
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it is possible to perform a further dimensional reduction of the phase space. First one 
notes that in vacuum and for K = the scale factor a only appears in eqs. (2.3)-(2.5) 
through the combination H = a/ a. One can then use the variables H and <fi and the 

Hamiltonian constraint to reduce the dimensionality of the phase space (^H, 0, <pj to 

two. It is appropriate to use the Hubble parameter because it is one of the cosmological 
observables. 

It is sometimes useful to rewrite eqs. (2.3)-(2.5) using the conformal time 77 defined 
by dt = a dr\ and the new variable 



a 1 



x =- = aH, (2.9) 



a 



where a prime denotes differentiation with respect to r\. In terms of x and r\ one has 



2 \<t> J \<t>) 2(2cu + 3)0 



• / , dV ^ T A , ,, S 2 duj 
*l + l$- 2V ) + M d*. 



87rp (^7 + 2) 

+K (2c. + 3)0a37-2- ^ n > 

In the rest of this paper we consider various specifications of the general scalar-tensor 
theory (1.2) corresponding to various combinations of the scalar field potential and of 
the cosmic fluid. 



3 Vacuum, free Brans-Dicke-like field, and any K 

In the case in which no material fluid is present and V = 0, eq. (2.10) can be written as 

uj (<f)') 2 — 6x<f)(f)' — 6<f) 2 [K + x 2 ) = . (3.1) 

Provided that uj 7^ 0, this is a second degree algebraic equation for 0' with reduced 
discriminant 

-j- = 3(f) 2 [(2uj + 3) x 2 + 2ujK] (3.2) 

and with solutions 

<f>'± (x, 0) = — ( x ± V 0) ) , (3.3) 
5 



^ 1 (x,0) = i{[2cu(0) + 3]a: 2 + 2^(0)} . 



(3.4) 



If uj > and K = or +1 it is always T\ > 0. If -K" = — 1 or uj < 0, there can be regions 
corresponding to T\ < which are forbidden to the dynamics. The use of eq. (2.11) 
yields 



x'± (x, 4>) = x 2 + ^- (x ± a/^i ) 



a; (2a; + 3) 



+ —[x± 



(x ± v^i") + # • (3.5) 



x' and 0' are determined once the values of x and are given. In general there are two 
values of both x' and 0' for a given pair (x, 0), corresponding to the upper and lower 
sign in eqs. (3.5) and (3.3), respectively. This property corresponds to the geometry of 
the phase space which is a two-dimensional surface embedded in the (x, x', 0, 0') space 
and consists of two sheets, each of which corresponds, respectively, to the upper or lower 
sign in eqs. (3.3) and (3.5) - we can call them "upper sheet" and "lower sheet". The 
(x, 0) plane is not the phase space but is a projection of it. As a consequence, the study 
of the (x, 0) plane can show projections of the orbits intersecting each other. Of course 
the true orbits living in the curved energy surface do not cross, due to the uniqueness 
of the solutions of the Cauchy problem for the system (2.3)-(2.5), but projections onto 
the (x, 0) plane of two different orbits located in the two different sheets can intersect. 

The orbits of the solutions can only change sheet (and they don't necessarily change 
in all the possible scalar-tensor scenarios) at points where the two sheets touch each 
other, i.e., where T\ = 0. This set of points forms the boundary of the forbidden region 



B = {(*,. 



T x (x, 0) = 0} 



(3.6) 



On points of B one has x' + = x'_ and 0' + = 0'_ = 3(f) x/uj. It should be noted that 
the phase space reduces to a plane instead of a curved two-surface when T\ vanishes 
identically, i.e., if K = and uj = —3/2. 

In the case of vacuum, V = 0, and K — the dimensional reduction proceeds by 
using comoving time and the variables H and - see Sec. 4 for an example. However 
if K = ±1 and comoving time is used, the scale factor cannot be eliminated from the 
field equations and the dimensional reduction cannot be achieved. To this purpose one 
must use conformal time when K = ±1. One notes that the potential can be neglected 
(V ~ 0) when the kinetic terms dominate the dynamics. In this case the evolution takes 
the universe away from the possibility of chaos and toward a two-dimensional phase 
space. If V((f)) dominates over the kinetic terms at late times, as in dark energy models, 
there is in principle the possibility of chaotic evolution, although this is not by all means 
guaranteed. 
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The fixed points of the system (3.3) and (3.5) correspond to (x', 0') = (0, 0); eqs. (3.3)- 
(3.5) then yield either = 0, which corresponds to infinite gravitational coupling and is 
unphysical, or x = and u> K = 0. Hence, 

• For K = the only fixed point is the Minkowski space corresponding to x — aH = 
0, which lies on the boundary B of the region forbidden to the dynamics. 

• For K — ±1 the only fixed points are of the form (x, 0) = (0, <j> ) where O is 
a root of the equation u (0) = 0. Such a fixed point corresponds again to a 
Minkowski space lying on B. If uo ^ everywhere there are no fixed points and, 
as a consequence, there are no limit cycles (a cycle must contain at least a fixed 
point). 



4 Vacuum, V = ra 2 </r/2, and K = 

In the situation of vacuum, a massive scalar 0, and K — one can use again the physical 
variables H and 0. Eq. (2.3) is rearranged as the algebraic equation for 



CO 
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(0) 2 -^00+(^0-^ 2 )0 2 = O (4.1) 



with discriminant 



and solutions 



A 2 (#,0) = ^- 



(4.2) 



ui m 2 



(2uj + 3)H 2 — 



± (tf,0) = ^(#0± v^) • (4.3) 

The phase space is again a curved two-dimensional surface embedded in the space 
(^H, 0, 0^ and composed of two sheets. The two sheets join on the boundary of the 
forbidden region described by A 2 = 0, which is equivalent to either = (unphysical) 



or 



^) = te J_M_. (4.4) 
;/ 6[2cj(0) + 3] V ; 



The fixed points (if o ,0o) are subject to the only restriction H = ±mA/0o/12 and are 
de Sitter spaces with constant scalar field O > located away from the boundary 
A 2 = 0. Again, if K ^ the scale factor cannot be eliminated and one cannot choose 
H as dynamical variable. For any value of K the wave equation for can be integrated 
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yielding either the trivial solution =constant (corresponding to general relativity with 
a cosmological constant) or 



J d(j) y/2u(<j>) + 3 = const. J ^ . (4.5) 

A complete description of the phase space is given by Santos and Gregory [10] for Brans- 
Dicke theory with oj =constant. They use the dynamical variables 



and conformal time r\. The only difference between the Brans-Dicke case and more 
general scalar-tensor theories are the terms in Co in the field equations of general scalar- 
tensor gravity. However, the dimensional reduction of the phase space cannot be achieved 
by using the variables X, Y, and rj. 

5 The case of vacuum, K = 0, and V ^ 

Let us consider again vacuum and a general self-interaction potential V(<f>) for the Brans- 
Dicke-like field 0, but let us impose the restriction that the universe has flat spatial 
sections (K = 0). This is the situation indicated by the Boomerang [23] and MAXIMA 
[24] experiments for our universe. In this section we use the physically transparent 
variables H and <fi and comoving instead of conformal time. The dynamical equations 
reduce to 



, dV nTr du ( -v 



2 \<f)J \<f>) 2(2cj + 3)0 

2cj + 3 ~d4>)^~ 2^~+ 3 V" ' r d(f) 



(5.2) 



+ f3tf + -JL_^ = -JL_( 2 y-^) . (5 . 3) 
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When uo ^ eq. (5.1) can be rearranged as the algebraic equation for 

uj(<py - 6#00 +(V- 6H 2 <f)) = (5.4) 

with reduced discriminant 

T 2 (H, 0) = [3 (2uj + 3) H 2 (f) - ujV] (5.5) 

and solutions 

± (if, 0) = -±- [3^0 ± ^2 1/2 (if, 0)1 • (5.6) 
u{(j>) L J 

The phase space is again a two-dimensional surface composed of two sheets corresponding 
to the lower and upper signs in eq. (5.6). Figures 1-9 illustrate the two-sheeted structure 
of the phase space for particular choices of c<j(0) and V((p). In general there is no 

guarantee that JF 2 > and there will be regions of the space ^H, 0, 0^ corresponding 

to JF 2 < forbidden to the orbits of the solutions. If K 7^ the reduction of the phase 
space to two dimensions cannot be performed because the scale factor appears in the 
expression K/a 2 in eqs. (2.3) and (2.4) and not in the form H — a/a. 

The phase space is flat if JF 2 = 0, i.e., if u — —3/2 and V = 0. In this case one has 

+ 2/70 = (5.7) 

and either is identically zero (a physically unacceptable solution) or oc a~ 2 . Then 
the Hamiltonian constraint is automatically satisfied and H = —H 2 , which yields the 
coasting universe with scale factor a = a (t — t ). 

In the general situation with 7^0 studied in this section, the equilibrium points 
correspond to (if, 0) = (H , O ) with H and 0o constants and must satisfy 




h = ±^-^- , o r ( ;-2r = o . (5.8) 

where V = V(0 O ) and Vq = dV/dcf)]^. 

6 The case 7 = 2/3 and V = 

Another special case corresponds to the choice 7 = 2/3 (or P = —p/3) and V = 0. In 
terms of the variables x and and of conformal time the field equations reduce to 

(9 +s f(9'-* + ^ 
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X 



x 2 - 



u(<j>) 





0' T , (00 du A 
_i_ 2r— A- K A — — 

^ ^ ^ 2 (2u + 3) <20 



where 



167rp (a; + 3) 



3 (2cj + 3) 

Eq. (6.1) can be rewritten as the algebraic equation for 0' 



7T (0') 2 - ^00' + 
6 



2 j. 2 



with discriminant 



A 3 (x, 0) 



(2uo + 3) x 2 - 2uo 



8vrp 



-X0 



and solutions 



/ ± (x,0) = -(a:0±v / A3 (^,0)) 



while eq. (6.2) yields 



x 2 - 



+ 



2cu0 2 
9 



± + ^ (x0 ± v^) + # 

do; 



2 

+ — ( £0 ± 

Co>0 



2lu 2 (2uj + 3) d0 



xd> ± 







(6.2) 
(6.3) 

(6.4) 

(6.5) 
(6.6) 



(6.7) 



The geometry of the phase space is again that of two sheets attached along the boundary 
A 3 (x, 0) = of a forbidden region and embedded in the space (x, 0, 0'). There are no 
fixed points in this case. In fact, the fixed points would correspond to constant Xq and 
0o and it is easy to see that eqs. (6.1) and (6.2) are incompatible when x and are both 
constant. 



7 Discussion and conclusions 

The geometry and dimensionality of the phase space of FLRW scalar-tensor cosmology, 
in which gravity is described by the action (1.2), can be quite complicated. In general 
one has two independent equations for scalar-tensor cosmology, the wave equation (2.5) 
for the scalar and eq. (2.4) for the scale factor a. Since these are second order equa- 
tions one is naturally lead to consider the phase space (a,d, 0,0). The Hamiltonian 
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constraint (2.3) provides a first integral and the orbits of the solutions lie on a hyper- 
surface in this space. Hence in general the phase space is a three-dimensional curved 

energy hypersurface embedded in the four-dimensional space (^a, a, 0, <j)j . When the spa- 
tial sections are flat (K = 0) the scale factor only enters the field equations through the 
Hubble parameter and it is convenient to choose as physical variables H and instead 
of a and 0. The Hamiltonian constraint (2.3) then forces the traiectories of the solutions 
to unfold in a two-dimensional hypersurface embedded in the three-dimensional space 

(jl, 0, <pj . However, there are also situations in which the universe has curved spatial 

sections (K = ±1) and one can again reduce the phase space to a two-dimensional sur- 
face embedded in a three-dimensional space. Table 1 summarizes the situations known 
in the literature and the results of the present paper. The latter are derived without 
specifying the coupling function uj{4>) and, in some cases, also for arbitrary scalar field 
potential V((p). The stationary points of the dynamical system are determined in this 
general situation. 

There is a simplification in the equations of scalar-tensor cosmology when the po- 
tential reduces to a mass term - previous works [9, 10] have taken advantage of this 
property. When a general scalar-tensor theory is considered instead of Brans-Dicke the- 
ory the time-dependence of u;(0) makes the dynamics more interesting even for vacuum, 
K = and a free scalar field. Even when the phase space is two-dimensional its geome- 
try is usually non-trivial. In general the phase space accessible to the trajectories of the 
solutions is composed of two sheets attached to each other along the boundary B of a 
region forbidden to the dynamics. The orbits of the solutions can change sheet only by 
passing through points of B. There are however scenarios in which the orbits are forced 
to stay in one sheet or the other but cannot change sheet - the sheet in which they live 
is decided by the initial conditions. The plane (H,<f)) (or (x, 0)) is only a projection 
of the true curved phase space and projections of the orbits can intersect in this plane 
(apart from the degenerate cases in which this plane reduces to the true phase space). 
In all the situations considered in this paper the fixed points of the dynamical system 
(when they exist) are found to be de Sitter spaces possibly degenerating into Minkowski 
spaces. Their stability is studied in Ref. [25]. 

A detailed discussion of the dynamics requires the specification of the functions 
a>(0) and V(4>). Although many forms of V{4>) have been considered in the literature 
on high energy physics and cosmology, there are very few suggestions about the form 
of the coupling function u>(4>), and they are usually dictated by purely mathematical 
considerations such as the ease of finding exact cosmo logical solutions (cf. Ref. [4]). The 
clarifications on the dimensionality and geometry of the phase space made here will be 
useful in detailed studies fixing the form of cj(0) and V((j>). 



11 



Acknowledgments 



It is a pleasure to thank Prof. Werner Israel for the hospitality at the University of 
Victoria, where this paper was completed. This work was supported by the Natural 
Sciences and Engineering Research Council of Canada (NSERC). 



References 

[1] C.H. Brans and R.H. Dicke, Phys. Rev. 124 (1961) 925. 

[2] RG. Bergmann, Int. J. Theor. Phys. 1 (1968) 25; R.V. Wagoner, Phys. Rev. D 1 (1970) 
3209; K. Nordvedt, Astrophys. J. 161 (1970) 1059. 

[3] C.G. Callan, D. Friedan, E.J. Martinec and M.J. Perry, Nucl. Phys. B 262 (1985) 593; 
E.S. Fradkin and A. A. Tseytlin, Nucl. Phys. B 261 (1985) 1. 

[4] V. Faraoni, Cosmology in Scalar-Tensor Gravity, Kluwer Academic, Dordrecht, 2004. 

[5] Y. Fujii and K. Maeda, The Scalar-Tensor Theory of Gravity, Cambridge University 
Press, Cambridge, 2003. 

[6] D. La and P.J. Steinhardt, Phys. Rev. Lett. 62 (1989) 376; E.J. Weinberg, Phys. Rev. D 
40 (1989) 3950; D. La, P.J. Steinhardt and E. Bertschinger, Phys. Lett. B 231 (1989) 
231; A.R. Liddle and D. Wands, Mon. Not. Roy. Astr. Soc. 253 (1991) 637. 

[7] F.S. Accetta and J.J. Trester, Phys. Rev. D 39 (1989) 2854; D. La, P.J. Steinhardt and 
E. Bertschinger, Phys. Lett. B 231 (1989) 231; R. Holman, EW. Kolb and Y. Wang, 
Phys. Rev. Lett. 65 (1990) 17; R. Holman, EW. Kolb, S. Vadas and Y. Wang, Phys. 
Lett. B 269 (1991) 252; A.M. Laycock and A.R. Liddle, Phys. Rev. D 49 (1994) 1827; J. 
McDonald, Phys. Rev. D 48 (1993) 2462. 

[8] T. Chiba, Phys. Rev. D 60 (1999) 083508; J.-P. Uzan, Phys. Rev. D 59 (1999) 123510; F. 
Perrotta, C. Baccigalupi and S. Matarrese, Phys. Rev. D 61 (2000) 023507; L. Amendola, 
Mon. Not. R. Astr. Soc. 312 (2000), 521; R. de Ritis, A. A. Marino and P. Scudellaro, 
Phys. Rev. D 62 (2000) 043506; X. Chen, R.J. Scherrer and G. Steigman, Phys. Rev. D 
63 (2001) 123504; G. Esposito-Farese, gr-qc/0011115; A. A. Sen and R.T. Seshadri, Int. 
J. Mod. Phys. D 12 (2000) 445; C. Baccigalupi, S. Matarrese and F. Perrotta, Phys. Rev. 
D 62 (2000) 123510; Y. Fujii, Gravit. Cosmol. 6 (2000) 107; Y. Fujii and T. Nishioka, 
Phys. Rev. D 62 (2000) 044011; O. Bertolami and J. Martins, Phys. Rev. D 61 (2000) 
064007; L. Amendola, Phys. Rev. D 62 (2000) 043511; Phys. Rev. Lett. 86 (2001) 196; 
B. Boisseau, G. Esposito-Farese, D. Polarski and A. A. Starobinski, Phys. Rev. Lett. 85 



12 



(2000) 2236; M.C. Bento, O. Bertolami and N.C. Santos, Phys. Rev. D 65 (2002) 067301; 
T. Chiba, Phys. Rev. D 64 (2001) 103503; A. Riazuelo and J.-P. Uzan, Phys. Rev. D 
62 (2000) 083506; G. Esposito-Farese and D.Polarski, Phys. Rev. D 63 (2001) 063504; 
N. Banerjee and D. Pavon, Class. Quant. Grav. 18 (2001) 593; V. Faraoni, Int. J. Mod. 
Phys. D 11 (2002) 471; D.F. Torres, Phys. Rev. D 66 (2002) 043522; F. Perrotta and 
C. Baccigalupi, astro-ph/0205217; A. Riazuelo and J.-P. Uzan, Phys. Rev. D 66 (2002) 
023535; S. Nojiri and S.D. Odintsov, Phys. Rev. D 68 (2003) 123512; hep-th/04050078; 

E. Elizalde et al, Phys. Lett. B 574 (2003) 1; E. Elizalde, S. Nojiri and S.D. Odintsov, 
hep-th/0405034; S. Nojiiri and S.D. Odintsov, Phys. Lett. B 595 (2004) 1; S. Matarrese, 
C. Baccigalupi and F. Perrotta, astro-ph/0403480; R. Catena et al., astro-ph/0406152; 

F. C. Carvalho and A. Saa, astro-ph/0408013. 

[9] S.J. Kolitch and D.M. Eardley, Ann. Phys. (NY) 241 (1995) 128; S.J. Kolitch, Ann. Phys. 
(NY) 246 (1996) 121. 

[10] C. Santos and R. Gregory, Ann. Phys. (NY) 258 (1997) 111. 

[11] D.J. Holden and D. Wands, Class. Quant. Grav. 15 (1998) 3271. 

[12] R.M. Wald, General Relativity, Chicago University Press, Chicago, 1984. 

[13] H. Nariai, Prog. Theor. Phys. 47 (1971) 1824; ibidem 48 (1972) 763; P. Chauvet and J.L. 
Cervantes-Cota, Phys. Rev. D 52 (1995) 3416; J.P. Mimoso and D. Wands, Phys. Rev. D 
51 (1995) 477; A. A. Coley, J. Ibanez and R.J. van der Hoogen, Int. J. Mod. Phys. D 38 

(1997) 5256; A. Billyard, A. A. Coley and J. Ibanez, Phys. Rev. D 59 (1999) 023507. 

[14] A.G. Riess et al, Astron. Astrophys. 116 (1998) 1009; S. Perlmutter et al, Nature 391 

(1998) 51; A.G. Riess et al., Astron. J. 118 (1999) 2668; S. Perlmutter et al, Astrophys. J. 
517 (1999) 565; A.G. Riess et al, Astrophys. J. 536 (2000); A.G. Riess et al, Astrophys. 
J. 560 (2001) 49; A.G. Riess et al., Astrophys. J. 607 (2004) 665. 

[15] A. A. Starobinsky, Gravit. Cosmol. 6 (2000) 157. 

[16] E. Gunzig et al, Mod. Phys. Lett. A 15 (2000) 1363. 

[17] V. Faraoni, work in progress. 

[18] S. Foster, Class. Quant. Grav. 15 (1998) 3485. 

[19] L. Amendola, M. Litterio and F. Occhionero, Int. J. Mod. Phys. A 5 (1990) 3861. 

[20] E. Gunzig et al., Class. Quant. Grav. 17 (2000) 1783. 

[21] E. Gunzig et al, Phys. Rev. D 63 (2001) 067301. 



13 



[22] A. Saa et al, Int. J. Theor. Phys. 40 (2001) 2295. 

[23] A.D. Miller et al, Astrophys. J. Lett. 524 (1999) LI; P. de Bernardis et al, Nature 400 
(2000) 955; A. E. Lange et al, Phys. Rev. D 63 (2001) 042001; A. Melchiorri et al, 
Astrophys. J. Lett. 536 (2000) L63. 

[24] S. Hanany et al, Astrophys. J. Lett. 545 (2000) L5. 

[25] V. Faraoni, Phys. Rev. D 70 (2004) 044037; Phys. Rev. D 69 (2004) 123520. 



14 



Table caption 



Table 1: Situations in which the geometry and dimension of the phase space are known 
(BD = Brans-Dicke theory, NMC = nonminimally coupled scalar field theory). 



Figure captions 

Figure 1: The upper sheet of the phase space of Brans-Dicke theory described by 
uj = 10 and scalar field potential V((f>) = m 2 (f) 2 /2 (for unit mass m). 

Figure 2: The lower sheet corresponding to the situation of fig. 1. 

Figure 3: The phase space of Brans-Dicke theory with massive scalar resulting from 
the two sheets of fig. 1 and fig. 2 glued together. The region forbidden to the dynamics 
appears as a hole that cannot be penetrated by the orbits of the solutions. 

Figure 4: The upper sheet of the phase space of the scalar-tensor theory described 
by uj = |0 — 1| _1 and scalar field potential V(4>) = m 2 2 /2 (for unit mass m). 

Figure 5: The lower sheet corresponding to the situation of fig. 4. 

Figure 6: The phase space resulting from the two sheets of fig. 4 and fig. 5 glued 
together. There is a forbidden region and the two sheets touch each other only at the 
boundary of this region. 

Figure 7: The upper sheet of the phase space of the scalar-tensor theory described 
by uj = (0 — l) -4 and scalar field potential V{4>) = V cos ((f) /m) (for m and V set to 
unity) . 

Figure 8: The lower sheet corresponding to the situation of fig. 4. 

Figure 9: The phase space resulting from the two sheets of fig. 4 and fig. 5 glued 
together. There is a forbidden region and the two sheets touch each other only at the 
boundary of this region. 
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